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ABSTRACT 

We investigate the possibility of detecting baryon acoustic oscillation features in the cosmic 
matter distribution by 3d weak lensing. Baryon oscillations are inaccessible even to weak 
lensing tomography because of wide line-of-sight weighting functions and require a special- 
ized approach via 3d shear estimates. We quantify the uncertainty of estimating the matter 
spectrum amplitude at the baryon oscillations wave vectors by a Fisher-matrix approach with 
a fixed cosmology and show in this way that future weak lensing surveys such as EUCLID 
and DES are able to pick up the first four wiggles, with EUCLID giving a better precision 
in the measurement. We also provide a detailed investigation of the correlation existing be- 
tween errors and of their scaling behavior with respect to survey parameters such as median 
redshift, error on redshift, error on the galaxy shape measurement, sky coverage, and finally 
with respect to the number of wiggles one is trying to determine. 
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1 INTRODUCTION 

Baryon acoustic oscillation (BAO) features are modulations in the 
cosmic matter distribution on very large spatial scales of roughly 
~ 100 Mpc//i (for a review, see Bassett & Hlozek 2010). These 
BAOs are the imprint of oscillations of the photon-baryon fluid 
in the early universe on the matter density field driven by grav- 
ity and the equation of state providing a restoring force, and they 
are observable in two primary channels: through the observation 
of anisotropies in the cosmic microwave background (CMB) and 
through galaxy surveys. The most important features such as their 
spatial scales, their signature in the CMB, their statistical prop- 
erties, and their dependence on cosmological parameters is very 
well understood analytically (Seljak & Zaldarriaga 1996; Hu & 
Sugiyama 1996; Montanari & Durrer 201 1 ; Sutherland 2012). 

Concerning the determination of cosmological parameters, it 
is a fortunate situation that they are observable at high redshifts 
through the primary CMB and at much lower redshifts in the galaxy 
distribution. Due to the fact that BAOs provide a standard yardstick 
at two different cosmological epochs, it is possible to constrain the 
density parameters of cosmic fluids and the possible time evolution 
of their equation of state parameters in a geometric way, breaking 
degeneracies that may arise if the equations of state of cosmological 
fluids are allowed to change with time. 

CMB observations earned out by the Cosmic Background Ex- 
plorer (COBE, Bennett et al. 1994; Wright et al. 1996) first re- 
vealed anisotropies in the CMB, but only the Wilkinson Microwave 
Anisotropy Probe (WMAP, Hinshaw et al. 2003, 2007; Nolta et al. 
2009; Larson et al. 20 1 1 ) had sufficient angular resolution such that 
the BAO scale of ~ 100 Mpc/A could be resolved at a comoving 
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distance of ~ 10 Gpc/h, revealing temperature modulations of the 
CMB of the order AT/T CMB 10" 5 at an angular scale of roughly 
2°, with subsequent higher harmonics. Likewise, galaxy surveys 
have now reached sufficient depth and solid angle that BAO fea- 
tures could be detected as modulations of the galaxy density of the 
order 10% in both radial and transverse directions. With the as- 
sumption of a galaxy biasing model, the longest wavelength BAO 
modes survive nonlinear structure formation to the present epoch 
(Meiksin et al. 1999) and will be targeted by future surveys for the 
precision determination of cosmological parameters (Dolney et al. 
2006; Angulo et al. 2008; Labatie et al. 2012), in particular dark 
energy (Seo & Eisenstein 2003; Eisenstein et al. 2007). Both av- 
enues have contributed significantly to the estimation of cosmolog- 
ical parameters and to the selection of most plausible cosmological 
models. 

Specifically, there are quite a number of detection reports with 
ongoing surveys, for instance with the Sloan Digital Sky Survey 
(SDSS, Eisenstein et al. 2005; Padmanabhan et al. 2007; Percival 
et al. 2007, 2010; Kazin et al. 2010; Padmanabhan et al. 2012; 
Mehta et al. 2012), the 2-Degree Field Galaxy Redshift Survey 
(2dFGRS, Percival et al. 2007; Beutler et al. 2011), the WiggleZ- 
survey (Parkinson et al. 2012) and Lyman-or data (Busca et al. 
2012) with subsequent determination of cosmological parameters 
which confirm spatial flatness and the low matter density found 
by CMB observations, if flatness is assumed prior to the analysis. 
Recent studies (Parejko et al. 2012; Zhao et al. 2012) were able 
to constrain neutrino masses. BAO modulations have been found 
as longitudinal as well as transverse modes in the galaxy density 
(Gaztanaga et al. 2009; Kazin et al. 2010; Gaztanaga et al. 2009) 
and their issues of model selection and parameter estimation have 
been addressed thoroughly (Cabre & Gaztanaga 201 1). 

The motivation for this paper is the fact that the detection of 
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BAOs as a modulation feature in the galaxy field depends on the 
assumption of a biasing mechanism (Gaztanaga et al. 2009; Des- 
jacques et al. 2010) which relates the galaxy number density to the 
ambient density of dark matter is well as a control of redshif t space 
distortions effects (Nishimichi et al. 2007; Taruya et al. 2009) and 
it would be desirable to measure the dark matter density directly. 
Weak lensing would be a prime candidate for such a measurement, 
but the wide line-of-sight weighting functions cause the weak lens- 
ing signal to depend rather on the integral of the spectrum of cold 
dark matter (CDM) than on individual, localized features, even in 
the case of tomographic lensing surveys (Hu 1999). This is the 
reason why investigate the sensitivity of 3d weak lensing (3dWL, 
Heavens 2003) for constraining the dark matter spectrum on BAO 
scales: 3dWL provides a direct estimate of the 3-dimensional mat- 
ter distribution and gives Gaussian errors on the amplitude of the 
CDM spectrum in wavelength bands from sparsely sampled data 
(Leonard et al. 2012). In this way, we aim to quantify the statisti- 
cal precision at which 3dWL constraints the CDM spectrum at the 
BAO wavelengths, and the statistical significance for inferring the 
presence of one or more wiggles from 3dWL data relative to the 
null-hypothesis of absent wiggles. 

After a short compilation of basic results concerning dis- 
tances, structure growth, structure statistics, and conventional weak 
lensing in Sect. 2, we recapitulate the main results of 3dWL in 
Sect. 3 and motivate its usage in constraining BAO wiggles. Our 
statistical approach and the estimation of statistical errors on the 
BAO measurement is given in Sect. 4, followed by a discussion of 
our main results in Sect. 5. 

The reference cosmological model used is a spatially flat 
vi'CDM cosmology with Gaussian adiabatic initial perturbations in 
the matter distribution. The specific parameter choices are Q m = 
0.25, n s = l,cr 8 = 0.8 and H = 100 h km/s/Mpc, with h = 0.72. 
The dark energy equation of state is set to w = -0.9 and we as- 
sume the dark energy to be smooth. The baryon density flj, = 0.04 
is used for correcting the CDM shape parameter and for predicting 
BAO-wiggle amplitudes and wave-vectors. 



2 COSMOLOGY AND STRUCTURE FORMATION 
2.1 Dark energy cosmologies 

In spatially flat dark energy cosmologies with the present matter 
density £l m , the Hubble function aH(a) = da/dt is given by 



(1) 



H 2 (a) CI,,, j f 1 

— p = -j- + (1 - fi„,)exp|3 J dlna(l+w(a)) 

with the dark energy equation of state w(a). A constant value w = 
- 1 corresponds to the cosmological constant. The relation between 
comoving distance x an d scale factor a is given by 
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with the Hubble distance xh = c/Hq as the cosmological distance 
scale. Redshift z and comoving distance are related by dz/dx = 
H(z)/c . 



2.2 CDM power spectrum 

The linear CDM-density power spectrum P(k) describes the fluctu- 
ation amplitude of the Gaussian homogeneous density field 6, 



with the the spectral index n s and the transfer function T{k). The 
restoring force provided by the baryon-photon fluid in the early 
Universe generates a set of wiggles in the spectrum P{k) and an 
overall suppression due to diffusion. Both effects are discussed in 
detail by Eisenstein & Hu (1998) and Eisenstein & Hu (1999) who 
also provide a fitting formula for T{k) in terms of the density pa- 
rameters £l„„ fl;,, and the Hubble parameter h. 

The spectrum P{k) is normalized in such a way that it exhibits 
the variance <r\ on the scale R = 8 Mpc//7, 
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with a Fourier transformed spherical top hat filter function, W(x) = 
3ji(x)/x, where j e (x) is the spherical Bessel function of the first 
kind of order I (Abramowitz & Stegun 1972). 

2.3 Structure growth 

The growth of density fluctuations in the cosmic matter distribu- 
tion can be described as a self-gravitating hydrodynamical phe- 
nomenon, in the limit of Newtonian gravity. Homogeneous growth 
of the density field 



Six, a) = D + (a)S(x,a = 1) 



(5) 



in the linear regime \S\ <k 1 is described by the growth function 
D + (a), which is the solution to the growth equation (Turner & 
White 1997; Wang & Steinhardt 1998; Linder & Jenkins 2003), 

d ' T D + (a) + -(3 + ■%+(«) = -^fi,„(a)D + (a). (6) 
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Nonlinear structure formation leads to a strongly enhanced 
structure growth on small scales, generates non-Gaussian features 
and, most importantly, wipes out BAO wiggles as features in the 
initial matter distribution. This can be understood in an intuitive 
way as corrections to the CDM spectrum in perturbation theory to 
order 11 assume the shape of integrals over polyspectra up to or- 
der 2n (which separate into a product of n spectra by application 
of the Wick's theorem, see the review by Bernardeau et al. 2002) 
and are therefore becoming insensitive to localized features that are 
not strongly influencing the normalization of P(k) (Springel et al. 
2005; Jeong & Komatsu 2006; Pietroni 2008; Matarrese & Pietroni 
2008; Crocce & Scoccimarro 2008; Nishimichi et al. 2009; Jeong 
& Komatsu 2009; Jurgens & Bartelmann 2012; Anselmi & Pietroni 
2012). Since nonlinear structure formation affects small scales first, 
we will target BAO wiggles with 3dWL beginning at the largest 
wavelength before proceeding to successively shorter wavelengths. 

2.4 Weak gravitational lensing 

The weak lensing convergence k provides a weighted line-of-sight 
average of the matter density 6 (for reviews, see Bartelmann & 
Schneider 2001; Munshi et al. 2008; Hoekstra & Jain 2008; Bartel- 
mann 2010), 
rXH 

k= d X W K ( X )S, (7) 
Jo 

with the weak lensing efficiency W K (x) as the weighting function, 
3£l m D + 



2 X h a 



(8) 



(S(k)S(k'y) = (2;r) 3 <5 D (A: - k')P(k) oc k"'T 2 (k), 



(3) 



and the lensing efficiency weighted galaxy redshift distribution, 
rewritten in terms of comoving distance. 
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G(x) 



(9) 



n(z) denotes a common parametrization of the redshift distribution 
of the lensed background galaxy sample, 



«(z) = ™o I — 



exp 



with l = 5> r (5 



(10) 



which can be rewritten in terms of a distribution in comoving dis- 
tance with the relation n(x)dx = n{z)dz using dx/dz = c/H(a). 
These expressions allow to carry out a Limber projection (Limber 
1954) of the weak lensing convergence, which yields the angular 
convergence spectrum C K (f), 
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(x)P(k = tlx). 
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We will formulate our derivations in terms of the lensing conver- 
gence k instead of the observable shear y because it is a scalar 
quantity and possesses identical statistical properties. Eqn. (11) il- 
lustrates why line-of-sight averaged weak lensing spectra are in- 
effective in picking up BAO wiggles (and is almost a repetition 
of the previous argument why nonlinear structure formation de- 
stroys BAO features): They provide only an integrated measure of 
the CDM spectrum P(k) weighted with wide weighting functions 
Wk(x) that is very insensitive to local features of the spectrum such 
as BAO wiggles. This argument holds even for advanced tomo- 
graphic surveys (Hu 1999; Takada & Jain 2004) and motivates the 
need of a 3-dimensional mapping of the cosmic matter distribution. 
With reference to Gaztafiaga et al. (2009) and Kazin et al. (2010), 
we would like to emphasize that weak lensing, due to its sensitiv- 
ity to gravitational shear components perpendicular to the line of 
sight, will provide measurements of BAO wiggles in the transverse 
direction. 



3 3D WEAK LENSING 

The method of 3dWL was introduced by Heavens (2003), who pro- 
posed to include distances of lensed galaxies estimated from their 
photometric redshifts to infer the 3-dimensional unprojected tidal 
shear, i.e. the second derivatives of the gravitational potential per- 
pendicular to the line-of-sight from distortions in the galaxies' el- 
lipticity. Therefore, this approach differs from estimations of the 
angular line-of-sight averaged spectrum C K (£) or corresponding to- 
mographic spectra C'J({) in the important respect that the statis- 
tics of the full 3-dimensional matter distribution is inferred without 
any averaging of shears with the line-of-sight galaxy distribution, 
which has been performed in eqn. (9). As such, 3dWL is particu- 
larly suited for the problem at hand, namely to provide a precise 
estimate of the amplitude of the dark matter power spectrum at 
the BAO wavelengths. Additionally, Heavens (2003) showed that 
if 3dWL is used for constraining P(k) at a fixed cosmology, the 
smallest errors are expected in the BAO regime of the CDM spec- 
trum. 

In this section, we recapitulate the main results of 3dWL in 
terms of the weak lensing convergence in the Fourier-convention 
we prefer to work with; please also refer to Castro et al. (2005), 
Massey et al. (2007), Heavens et al. (2006), Kitching et al. (2008) 
for a detailed description of the theory, to Munshi et al. (201 1) for 
higher-order statistics through 3dWL and to Ayaita et al. (2012) for 
details of our numerical implementation. We assume spatial flat- 
ness and lensing in linearly evolving structures, which can be, in 



principle, relaxed from the 3dWL point of view (Pratten & Mun- 
shi 2013). The impact of systematic errors is nicely investigated by 
Kitching et al. (2008), and for an application to observational data 
we refer the reader to Kitching et al. (2007). 

The most natural choice for carrying out a Fourier transform 
in spherical coordinates is a combination of spherical harmonics 
for the angular and spherical Bessel functions for the radial depen- 
dence. We can therefore write the transformation for the conver- 
gence k as 



(see Ballinger et al. 1995; Heavens & Taylor 1995), where je and 
Y fm are, respectively, a spherical Bessel function of the first kind 
and a spherical harmonic, and = (6, (f). There exist algorithms for 
fast computation of Kc,„{k) (Percival et al. 2004; Rassat & Refregier 
2012; Lanusse et al. 2012; Leistedt et al. 2012). Such a transforma- 
tion is particularly convenient as the combination of j( and Y( m is 
an eigenfunction of the Laplacian in spherical coordinates, leading 
to a quite simple relationship between the coefficients of the density 
field 6( m (k) and the lensing convergence K[„,(k) as the observable: 
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with implicit assumption of the Einstein summation convention 



X(k,k')Y(k',k") : 
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k"dk' X(k,k')Y(k',k"). 



(15) 



It is then possible to construct an estimator for x tm (k) by in- 
cluding the uncertainty of the galaxy distance estimates coming 
from errors in the measurements of redshift. If we denote by x the 
true radial coordinate of a galaxy, and by x' the one inferred by its 
observed redshift z' = z(x'), then they will be related by the proba- 
bility p(x'W)> which we assume to be Gaussian for simplicity: 



p(x'k)dx : 



1 



: exp 



(z(x)~z(x')) 2 



2cr? 



dz' 



(16) 



where <x- is the width of the distribution and is assumed to be con- 
stant throughout the entire galaxy sample. Furthermore, galaxies 
receive a statistical weight according to their distribution in dis- 
tance n(x)dx- Following the derivation in Heavens (2003), we de- 
fine the two additional matrices 

Z e (kX) = lfx' 2 dx' fdxp(x'\x)jdk'x)jf(kx'l (17) 
M((k, k') = 2 Jx 2 dxn(x)je(kx)j f (k'x), (18) 

where nix) ls the number density of galaxies, as defined in 
eqn. (10). These matrices describe the correlations in spheri- 
cal Fourier modes generated by the measurement process: While 
rj[(k,k') describes mode couplings due to weak lensing, Z[(k,k') 
and M[(k,k') define, respectively, the contributions in the mode 
couplings coming from redshift errors and from the galaxy distri- 
bution along the radial coordinate x- 

We restrict ourselves to observations of the entire sky. In this 
case, the expression for the estimator k [m of the convergence is then 
expected to be 
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_ 3£l m W+l) B ( (k,k") 



(19) 



where the mode-coupling matrix B t (k,k') describes two integra- 
tions over k\ and ki\ 



B f (k, k") = Z e {k,ki)Me(kuk 2 )ne{k2X)- 



(20) 



Since the average values of a field like K{,„(k) are zero for all- 
sky surveys, we can only infer information about any parameter the 
field may depend on by means of its covariance, 



(KC,n(k)k* m (k')) = S K ,f(k,k') + N K ,f(k,k') = C^Kk') 



(21) 



which consists of a signal term S K ((k,k') and a noise term 
N K f(k,k'). The signal term S K j can be calculated directly from 
eqn. (19): 



3Q„ 



t(f + 1) 



with the abbreviations 



B t ik,k") B ( (k',k") 



ik"f 



(k"Y 



B t {k,k") = Z e ik,h)M ( ikuk 2 )iliihX) 
B c (k',k") = Z ( (k',k 3 )M ( (k3,k 4 )r lc (k 4 ,k") 



Pdk"), (22) 



(23) 
(24) 



with implicit integration over k t , k 2 and k 3 , k 4 . The corresponding 
noise part N Kt ( is given by 



#„(*.*') = -fMdKk'), 



(25) 



which is proportional to the shape noise a 2 e , namely the variance of 
the galaxy ellipticity distribution. It is important to notice that N K [ 
is independent of cosmology or variations in the CDM spectrum 
P(k). Intrinsic ellipticity correlations were neglected, which would 
greatly complicate the 3dWL description. 




k [h Mpc 



Figure 1. Ratio between the power spectrum with baryon acoustic oscil- 
lations P(k) and the smooth power spectrum P s (k). The largest wiggles in 
amplitude have been highlighted and labeled as w a ,a = 1,..., 6. These are 
the wiggles used to parametrize the power spectrum in our Fisher matrix 
approach. 



Fisher matrix specifies what are the best errors to expect for the 
inferred parameters 8 a via the Cramer-Rao relation. 

It can be proved that, since f -measurements are independent 
in the case of full-sky coverage, we can reformulate eqn. (28): we 
consider our estimator to be k{„,(k) and its covariance as defined in 
eqn. (21), and find 



r (max 

Fcfi = -j- 2 Q-t + l)tr[(C,/-' d a C K , t ) x (C, 



/ 1 3/jCV/) 



(30) 



4 DETECTING BAO WIGGLES 
4.1 Construction of the Fisher matrix 

We choose a Fisher matrix approach to determine how precisely 
3dWL can constrain baryon acoustic oscillations in the matter 
power spectrum P(k). The Fisher matrix is a square matrix whose 
elements are defined as the expectation values of the second deriva- 
tive of the logarithmic likelihood with respect to the fiducial param- 
eters 6 a and Qp\ 



a/3 



I d 2 \n£ 

\dO a dG p 



(26) 



As a general statement, if the likelihood £, can be expressed as an 
A'-dimensional Gaussian 



£ = 



exp — x C x 

^(2tt) n det(C) \ 2 



(27) 



where x is a generic data vector and C is the corresponding covari- 
ance, we can then write 



F aj9 = Itr[(C- 1 fl a C)x(C- I ^C)], 
or, equivalently, 

F a fs = -tr^lnCx^lnC], 



(28) 



(29) 



where d a and dp stand for the derivatives with respect to the pa- 
rameters 6 a and 8p. Given a particular experimental framework, the 



In our specific case, we consider the power spectrum as 
parametrized not by usual cosmological parameters, such as Q.,„ or 
<r 8 , but rather by its own wiggles amplitudes, namely by the values 
assumed by P(k) in a range of k where the wiggles w a are located. 
In Fig. 1 we plot the wiggle-only power spectrum, i.e. the ratio be- 
tween the power spectrum with baryon acoustic oscillations P(k) 
and an equivalent, smoothed out spectrum that has the same shape 
as P{k) but shows no oscillating feature, P s (k). We highlight the 
wiggles that have been used to parametrize P(k), w„,a = 1, 6. 
By calculating the derivative d a C of the covariance with respect to a 
variation of the amplitude of a maximum number of wiggles « w , we 
can build up (eqn. 28) a Fisher matrix F„p, where a,B = 1, ...,n w . 
Such a matrix carries information about the best errors to expect 
on the detection of each wiggle w a , a = 1, ...,h w , and the cross- 
correlations between inferred wiggle amplitudes. Given our aim, 
what we are actually performing in the calculation of d a C is a func- 
tional derivative, also known as Frechet derivative. In fact, we can 
imagine the power spectrum as depending on features, i.e. the wig- 
gles in Fig. 1. Each one of them can be approximated to a sin-like 
function defined in a range of k as wide as A/2, where A is the 
wavelength of the function itself. The covariance derivative is nu- 
merically estimated for one wiggle at a time as a finite difference: 



d a C K s 



2e ' 



(3D 



where C*. are the covariance matrices calculated using the power 
spectra P^(k) and e is an arbitrarily small number. The spectra 
P%ik) are equivalent to the original P(k) for all k of the domain, 
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Figure 2. The ratio between P(k) and P s (k). The figure shows a pictorial 
representation of the variation of the amplitude of one wiggle when calcu- 
lating the derivative of the covariance matrix, d a= 2Cg. 
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Table 1. Basic survey characteristics used for the Fisher-analysis: median 
redshift Zmed of the galaxy sample, galaxy density per squared arcminute ii, 
sky coverage fraction / s j; y , redshift error tr z and shape measurement error 
<t £ of the surveys EUCLID, DES and a hypothetical deep-reaching survey 
labeled DEEP. 

exception made for the wave numbers belonging to the interval I a 
that corresponds to wiggle w„. In this interval P^(k) is then 

P%k) = P(k) ± eP(k). (32) 

We would like to point out that, since what we are actually per- 
forming by means of the spectrum variation in eqn. (32) is in a way 
a logarithmic derivative of Ci, the denominator in eqn. (31) lacks a 
factor P{k) and is therefore just two times the fraction of the spec- 
trum used in the variation. In Fig. 2 we show a representation of an 
example of the variation performed in the calculation of the deriva- 
tive of the covariance matrix, d a C; in this case the second wiggle, 
W2, has been considered. 

It is worth noticing that the Fisher-matrix approach for infer- 
ring the error <r a on the dark matter spectrum P(k a ) (where k a are 
simply the k e /„, for brevity) as a Gaussian standard deviation 
is perfectly justified because of the linearity of the lensing observ- 
able and the linearity of the random field, so we do not need to use 
Monte-Carlo sampling for evaluating the likelihood £.(P(k a )) and 
to measure its widths <r a from Monte-Carlo samples of the likeli- 
hood. 

As noise sources for the inference of P(k a ), we consider a 
Gaussian shape measurement error cr e for the galaxy ellipticities, 
which are assumed to be intrinsically uncorrelated, and a Gaus- 
sian error a. for the redshift determination uncertainty. Likewise, 
we work in the approximation of neglecting all geodesic effects 
(Seitz & Schneider 1994; Seitz et al. 1994) like deviations form the 
Born approximation, lens-lens couplings (Shapiro & Cooray 2006; 
Krause & Hirata 2010), source clustering (Schneider et al. 2002), 



Figure 3. The marginalized errors on the first four wiggles, normalized with 
respect to their amplitude, as the maximum t used for the calculation of 
the Fisher matrix increases. It shows that choosing ( mllK > 100 brings no 
particular advantage in the precision of the wiggles measurement. In the 
calculation of the errors we considered the power spectrum as parametrized 
by the first 4 wiggles. 

source-lens correlations (Hamana et al. 2002), and deviations from 
Newtonian gravity (Acquaviva et al. 2004). While performing the 
necessary variations for computing the Fisher matrix, we keep all 
other cosmological parameters fixed and calculate everything us- 
ing an f -range between £ m ; n = 2 and f max = 100 (please see the 
next section for a justification of this choice). As surveys, we con- 
sider the cases of EUCLID, DES, and a hypothetical deep-reaching 
galaxy survey we refer to by the name DEEP. The respective survey 
properties are summarized in Table 1 . All the details concerning the 
numerical implementation used in the calculation of the covariance 
matrices can be found in Ayaita et al. (2012). 

4.2 Statistical errors 

The error cr a for inferring the amplitude of the CDM spectrum 
P(k a ) at wiggle positions k a is given by the Cramer-Rao relation, 

crl = (F-'U, (33) 
and 

oi = 1/F„, (34) 

for marginalized and conditional likelihoods, respectively. 

Before carrying out our analysis for surveys like EUCLID, 
DES, and DEEP, we implement some tests in order to determine 
the optimal value for the maximum number of modes to be used in 
the calculation of the Fisher matrix, C m . dx . Besides, we tried to find 
out how and by how much are the errors sensitive to some of the 
usual survey parameters, such as 

(i) the shape noise <r € ; 

(ii) the error or. in the measurement of redshift; 

(iii) the median redshift z mcd ; 

(iv) the fraction of sky coverage f S ty. 

Throughout these tests, when not stated otherwise, we make use of 
a default set of survey parameters such that <x- = 0.02, cr E = 0.3, 
Zmcd = 0.9, / s k y = 0.4 and h = 20. Additionally, we assumed we 
want to constrain simultaneously the first four wiggles (see Fig. 1). 
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Figure 4. The relative marginalized error for the first 4 wiggles, as the shape 
noise cr E varies. Again, we assumed tr- = 0.02, <r E = 0.3, z m ed = 0.9, 
fskr = 0.4 and n = 20. As we expected, larger values of cr E bring along 
larger errors on the detection of the wiggles. 



Figure 5. The relative marginalized error for the first 4 wiggles as a func- 
tion of the eiTor in the measurement of the photometric redshift, cr z . The 
sensitivity of the errors on the value of o% is not as steep as it is for the 
shape noise cr E . In addition, the increasing rate of cr a seems to change as 
we consider higher order wiggles. 



We start our investigation by determining the errors cr a , a = 
1, 2, 3, 4, as the maximum number of modes C mla in the summation 
in eqn. (30) increases. Please refer to Fig. 3 for a plot of the be- 
havior of <r„ normalized to the oscillation amplitude A„, where A a 
is defined as the maximum value of \P(k) - P s (k)\ for each wiggle. 
In particular, we considered ( m!ix = 10, 30, 100, 300, and observe 
that, after C mia = 100, there is practically no gain in the precision 
with which the first 4 wiggles would be constrained in a 3dWL 
approach. Therefore, we decide to stick to a maximum number of 
modes of 100 for all the subsequent calculations. This is a fair ap- 
proximation also from a theoretical point of view: in fact, extending 
too much the £ interval for the F a p summation could make us fall 
out of the linear regime; in addition, the assumption of a Gaussian 
shape for the likelihood £, could not to be anymore reasonable in 
such a multipole range (Heavens 2003). 

In Fig. 4, we show what happens as soon as we keep all the 
survey parameters fixed and vary the shape noise cr e . As one can 
expect, larger values of cr e lead to larger <r a , therefore to a greater 
uncertainty in the detection of the wiggles. Moreover, the rate at 
which cr a grows with the shape noise seems to be of the type of 
a power law, and seems independent of the wiggle considered, at 
least for the wiggles sample we evaluated. The situation turns out 
to be similar when the error in the determination of redshift o~- is 
considered, in Fig. 5: increasing a- still produces larger errors on 
all the wiggles under investigation, although here the relation is 
somewhat slower, especially as long as <r z < 0. 1 ; the relation also 
appears to be slightly dependent on the wiggle, becoming steeper 
as higher order oscillations are taken. In fact, by incrementing the 
error on redshift from 0.01 to 0.1, we get an error larger only by 
a factor of ~ 2 on the first wiggle and by a factor of ~ 8 on the 
second. Additionally, the correspondence between a a and the me- 
dian redshift of the survey (Fig. 6) seems again like a power law 
that gives larger errors for an increasing Zme&, and is independent of 
the wiggle. Qualitatively, this trend makes sense in light of the fact 
that, as we increase the median redshift, we keep fixed all other sur- 
vey parameters such as, for example, the galaxy density per squared 
arcminute h. By doing so, we consider surveys where a number h 
of galaxies is distributed over a deeper cone, meaning that we are 
actually sampling the 3D convergence field in a more diluted way, 
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Figure 6. Relative marginalized errors on the first 4 wiggles when the me- 
dian redshift z. m <:d of the survey increases. Since all other parameters are 
kept fixed, especially the mean galaxy density per squared arcminute, h, in- 
creasing Zmed is equivalent to having more and more diluted surveys, where 
the same amount of galaxies is distributed along a deeper survey cone. 



and therefore inheriting a larger noise. Naturally, varying the sky 
covarage propagates to the errors cr„ cc 1 / Jf^y. 



4.3 Detectability of BAO wiggles 

In this section, we would like to present the results obtained when 
estimating the best errors to expect on BAO wiggles for the surveys 
EUCLID, DES, and DEEP (please see Table 1 for specifications). 

We started our analysis by calculating both marginalized 
(o~ a = -\/(F~ l )„„) and conditional errors (cr„ =1/ VF„„) relative 
to the wiggle amplitude. These errors were computed for the three 
types of surveys, considering the first four oscillations, as shown in 
Fig. 7. As we could expect, marginalized errors are always larger 
than the correspondent conditional ones, namely the <r a on each 
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Figure 7. Conditional (magenta lines) and marginalized relative errors 
(black lines) for the three surveys under investigations EUCLID (solid line), 
DES (dashed line) and DEEP (dash-dot line) as a function of the wiggles, 
when the first 4 oscillations are simultaneously constrained. 



Figure 9. Fisher confidence ellipses for DES, when one tries to simulta- 
neously constrain the first four wiggles. Again, contours areas correspond 
to 1 - 2 - 3cr and the number in every panel is the correlation coefficient; 
the axes represent variation of wiggles in terms of their amplitude fraction. 
Also in this case we can observe correlation between the amplitudes of P(k) 
at different wiggles positions. 
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Figure 8. Confidence ellipses for the first four wiggles in a EUCLID-like 
survey, showing that the wiggles are indeed highly correlated. In fact, tak- 
ing, for example, two consecutive wiggles, such as w% and W3, we see that 
by increasing the amplitude of P(k) at the position of wiggle 2, we must 
then have the amplitude at W3 decreased in order to remain in the confi- 
dence region. The x- and y-axes show the variation of the wiggle in terms 
of percentage of its amplitude A a , the three contours areas correspond to 
1 - 2 - 3<x, and every panel shows the correlation coefficient in the upper- 
right corner. 



wiggle when we assume to know precisely all the other wiggle am- 
plitudes. 

We continue the investigation considering the confidence el- 
lipses calculated from the corresponding Fisher matrices obtained 
for the three surveys. We assume we are aiming to jointly constrain 
the first four wiggles and plot the results in Figs. 8, 9, and 10 for, 
respectively, EUCLID, DES, and DEEP. The sizes of the ellipses, 
whose contours stand for 1 - 2 - 3cr, already tell us that, among the 
ones evaluated, EUCLID will probably be the survey with largest 
constraining power on the BAO wiggles. It is of particular interest 
noticing the orientation of the ellipses, or their correlation coef- 
ficients (upper-right corner in every panel), that tell us something 
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Figure 10. Confidence ellipses for the DEEP survey, with the variation of 
the oscillations in terms of the wiggle amplitude. Again, we assumed we 
wanted to constrain jointly the first four wiggles, contour areas stand for 
1 - 2 - 3<r and the correlation coefficient can be read in the panels. 



about the interdepence between different wiggles: in fact, neighbor- 
ing wiggles are anti-correlated , i.e. increasing the amplitude of the 
power spectrum in correspondence to one oscillation would cause 
the P{k) at the position of the adjacent wiggle to take smaller val- 
ues in order to stay among the confidence region, and vice versa, 
whereas the opposite holds for alternated wiggles. 

In order to better understand whether the constraining power 
of the three surveys will allow us to detect any oscillatory feature 
in the CDM power spectrum, we plot the a„ obtained from the 
Cramer-Rao relation in eqn. (33) as error bars in the usual wiggle- 
only power spectrum for DEEP (Fig. 13), DES (Fig. 12), and EU- 
CLID (Fig. 11). Since what is shown is a ratio between P(k) and 
a smooth spectrum, the cr a have of course also been normalized 
with respect to P s (fc). The four different panels show how the errors 
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Figure 11. The bars show the marginalized errors in the detection of the 
wiggles, normalized with respect to P s (k), when the first 2 (upper left 
panel), 3 (upper right panel), 4 (bottom left panel), 5 (bottom right panel) 
wiggles are used to parametrize the power spectrum. Here, we considered a 
EUCLID-like survey with cr : = 0.1, cr e = 0.3, ; med = 0.9, / sky = 0.5 and 
n = 30. 
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Figure 12. Marginalized errors on the detection of BAO wiggles when one 
tries to detect the first 2 (upper left panel), 3 (upper right panel), 4 (bottom 
left panel), 5 (bottom right panel) wiggles at the same time, for a DES-like 
survey (<t z = 0.02, tr e = 0.3, z me d — L5, /sky = 012 and n = 10). Again, 
we considered the relative marginalized errors, dividing by P s (k). 



Figure 13. The marginalized errors on the detection of the wiggles in the 
power spectrum, relative to P s (k), for the hypothetical survey DEEP, char- 
acterized by the parameters cr- = 0.02, tr e = 0.3, z mcd = 1.5, f&y = 0.1 and 
n = 40. 



EUCLID give a better performance, allowing for up to 4 wiggles to 
be simultaneously constrained, with EUCLID giving smaller errors 
overall (Fig. 11). 

A better comparison between the three surveys can be car- 
ried out analyzing Figs. 15 and 14, where we plotted relative errors 
cr a /A a as functions of the maximum number of wiggles we want to 
jointly constrain, n w , and we collate results coming from, respec- 
tively, DEEP and EUCLID, and DES, and EUCLID. It becomes 
straightforward that a DEEP-like survey cannot compete against 
EUCLID: the relative errors coming from DEEP are always larger 
than the latter's, independently of the maximum number of wig- 
gles to be constrained, and remain safely under the unity only for 
the first 2 wiggles. 

The situation is better for DES, which gives results quite sim- 
ilar to EUCLID's up to /7 W = 4, and an even higher performance 
afterwards, since its o~ a are increasing at a smaller pace with re- 
spect to EUCLID's. This is not very useful, anyway, as the relative 
errors become larger than 1 for « w > 4 for both surveys. 

Concluding, EUCLID seems to grant the best results, allowing 
for the simultaneous detection of up to 4 wiggles with expected 
errors that are smaller than the ones predicted for both DES and, of 
course, DEEP. 



change when we try to jointly constrain the first 2, 3, 4, or 5 wiggles 
with a 3dWL approach. 

What these and the following plots show, first of all, is an ex- 
pected feature: as we increment the number of wiggles we expect to 
simultaneously examine, the precision with which the amplitudes 
P(k a ) would be measured gets poorer and poorer for all the oscil- 
lations. Our purpose would then be to evaluate how many BAO 
wiggles one is allowed to constrain before the errors on them be- 
come too large. It can be seen that all three surveys would allow 
for quite good constraints on the first 2 wiggles. The hypothetical 
survey DEEP already shows error bars of the order of the wiggle 
amplitude A a when the first 3 wiggles are considered, and the er- 
rors become much larger than A a (a > 1) as soon as one tries to 
detect 4 or more wiggles (Fig. 13). On the other hand, DES and 



5 SUMMARY AND CONCLUSIONS 

Subject of this paper has been a statistical investigation on whether 
future weak lensing surveys are able to detect baryon acoustic os- 
cillations in the cosmic matter distribution by application of the 3d 
weak lensing method. For a fixed vvCDM cosmology, we have es- 
timated the statistical precision cr„ on the amplitude of the CDM 
spectrum P(k) at the BAO wiggle positions in a Fisher-matrix ap- 
proach. Throughout, we worked under the assumption of Gaus- 
sian statistics, independent Fourier modes and in the limit of weak 
lensing. Noise sources were idealized and consisted in indepen- 
dent Gaussian-distributed shape-noise measurements for the lensed 
background galaxy sample, as well as a Gaussian error for the red- 
shift determination. As surveys, we considered the cases of EU- 
CLID, DES and a hypothetical deep-reaching survey DEEP. 
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Figure 14. Relative marginalized errors (with respect to the wiggles am- 
plitude) as a function of the maximum number of wiggles that we try to 
simultaneously detect, n w . Here, we show the results for a EUCLID-like 
survey (solid lines) and a DES-like survey (dashed lines). Different colors 
and point types correspond to the different wiggles on which the error tr a 
is calculated. The relative errors coming from a DES-like survey turn out 
to be larger than EUCLID'S as long as we try to observe a maximum of 
four wiggles. DES gives a better performance for n w > 4, but the <r a on the 
measurement of the higher order wiggles prove to be too large with respect 
to their amplitude A a , anyway. 
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Figure 15. Relative marginalized errors as a function of n w for a EUCLID- 
like survey (solid lines) and a DEEP-like survey (dashed lines). The relative 
cr a for a DEEP survey appear to be always larger than the ones obtained 
from EUCLID, independently of the maximum number of wiggles n w we 
want to constrain. 

(i) We have constructed the Fisher matrix considering our model 
as parametrized by the amplitudes of the CDM power spectrum 
at the baryon acoustic oscillations anticipated positions. In partic- 
ular, we started taking the two BAO wiggles with largest ampli- 
tude and progressively increased the number of oscillations consid- 
ered. Keeping the cosmology fixed to a standard wCDM param- 
eter choice, we carried out variations of P(k) that preserved its 
wiggle-shape in those wave number intervals; we then estimated 
the Fisher matrix accordingly, in order to quantify whether the 
statistical power of future weak lensing surveys suffices to place 
bounds on the amplitudes of the considered harmonics. By means 



of the Cramer-Rao relation, we calculated the best errors a a to ex- 
pect for the amplitudes of P(k) at wiggle positions. 

(ii) The sensitivity of <x„ with respect to some typical survey- 
parameters was tested. In particular, we considered the shape noise 
a € , the redshift error cr ; , the median redshift z me d, and the sky cov- 
erage / sky . We found that, as expected, increasing the uncertainty 
in the estimate of either the redshift or the galaxy shapes brings a 
larger error in the inference of the presence of wiggles, and that the 
sensitivity of these errors on <t- is less pronounced for small values 
of cr ; , although it grows as soon as we consider higher order wig- 
gles or large <r z . An increase of z me d leads as well to larger errors 
on the wiggles amplitudes, as one would expect from considering 
less and less populated surveys; for the same reason, a wider sky 
coverage, i.e. larger f&y, yields to higher precision in constraining 
the wiggle amplitudes cr a oc 1 / -Jf^. 

Overall, we may conclude that the volume of a survey seem to be 
overcoming the importance of a high precision in the redshift mea- 
surement of galaxies, at least for cr z < 0.1 - 0.2. 

(iii) Finally, we evaluated the cr„ for the surveys under investi- 
gation and found that, among them, EUCLID gave the best results, 
potentially allowing for the detection of up to the first four BAO 
wiggles with a good statistical confidence. Given our tests on the 
sensitivity of the errors on P(k a ) to certain survey parameters, we 
may conclude that EUCLID'S good performance is probably due to 
the volume of the survey in terms of total galaxy number and sky 
coverage, that seem to prevail over the negative effects brought by 
the error on redshift measurements, cr_, quite larger than the ones 
predicted for the other two surveys. 

Given these results, we conclude that measurements of BAO 
wiggles based on future weak lensing data are entirely possible, and 
avoid issues related to galaxy biasing and redshift-space distortions. 
We forecast a detection of the first four wiggles with EUCLID and 
DES by applying 3dWL techniques. Future developments from our 
side include estimates of the precision that can be reached on in- 
ferring dark energy density and equation of state by including the 
estimate of the BAO scale at low redshifts probed by lensing to the 
estimates at intermediate redshift provided by galaxy surveys and 
those at high redshifts such as the CMB. Additionally, we are inves- 
tigating the impact of systematical errors on the estimation process 
from 3dWL-data and biases in the estimation of BAO-wiggle am- 
plitudes. 
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